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1. Introduction 

The crucial tool in the search for the quark-gluon plasma (QGP) is a probe 
to test if the strongly interacting medium produced in nuclear collisions consists 
of confined or deconfined quarks and gluons. In this survey, we will show how 
quarkonia can be used as such a tool. First of all, this requires an understanding 
of the dynamics of quarkonium production in the absence of a thermal medium, 
i.e., in hadron-hadron collisions. The theory of quarkonium production can today 
be tested on a variety of different states (J/ip, % c , ip', T, T', T"); it is confirmed 
by data over a vast range of collision energies, up to 1.8 TeV. Given quarkonium 
production dynamics, we then have to address three distinct problems to obtain 
a viable probe. What are the effects of confined and of deconfined matter on the 
production of the different quarkonium states? How can pre-equilibrium phenomena 
affect the production? We begin with a short summary of the answers which we 
shall obtain to these questions. 

1 . 1 Preview 

To be specific, we consider charmonium states; but the arguments are in gen- 
eral applicable to bottonium states as well. The theory of charmonium-hadron 
interactions predicts that the tightly bound charmonium ground state J/ip cannot 
be dissociated by hadronic matter at temperatures below 0.5 GeV. The resulting 
transparency of confined media to J/i/Vs can be tested experimentally in nuclear 
matter, by studying //^-production with a P6-beam incident on a hydrogen or 
deuterium target. In contrast, a quark-gluon plasma contains deconfined gluons, 
and these are hard enough to break up a physical J/ip. 

Hadron-nucleus collisions also provide the information needed to eliminate pos- 
sible initial state effects. By comparing different hard processes with and without 
possible final state interactions, the role of the nuclear medium on the initial state 
can be clarified and the final state effects can be singled out. 

Combining hadron-nucleus and nucleus-nucleus studies to determine the effects 
of confined and deconfined systems on quarkonium production, we can thus test 
colour deconfinement. In a nutshell: once initial state effects are removed, confined 
systems have no effect on the J/ip, but suppress the ip'; once deconfinement sets 
in, both J/ip and ip'aie suppressed, but the ip' more strongly. 

1.2 Probing Colour Deconfinement 

After this short preview of things to come, we return to the general problem of 
establishing quark-gluon plasma formation. In high energy nuclear collisions, two 
beams of partons collide; the partons are initially confined to the colliding nucleons. 
This confinement can be checked, e.g., by studying primary high mass Drell-Yan 
dilepton production and observing, except for possible nuclear shadowing effects, 



the same parton distribution functions as in deep inelastic lepton-nucleon collisions. 
After the primary collision, we expect abundant multiple interactions, leading to a 
rapid increase of entropy, quick thermalisation and hence the production of strongly 
interacting matter. The fundamental question is whether confinement survives this 
thermalisation. If it does, we have hadronic matter - if not, a quark-gluon plasma. 
We expect confinement to be lost if the parton density sufficiently surpasses that 
present in a hadron-hadron interaction, so that partons can no longer be assigned 
to specific hadrons. How can we check if this has happened? 

The QGP is a dense system of deconfined quarks and gluons. Its density is 
in fact the reason for deconfinement: in a sufficiently dense medium, the long- 
range confining forces become screened, so that only short-range (<g Aq^ d ) in- 
teractions between quarks and gluons remain. To study such a medium and de- 
termine its nature, we therefore need probes which are hard enough to resolve 
the short sub-hadronic scales and which can distinguish between confined and de- 
confined quarks and gluons. In addition, the probe must survive the subsequent 
evolution of the medium; therefore it certainly cannot be in equilibrium with the 
later stages of matter. Two hard, strongly interacting signals produced before 
equilibration and distinct from the medium have been proposed as probes for con- 
finement / deconfinement: 

- heavy quark-antiquark resonances (charmonium, bottonium) [1,2], and 

- hard quarks or gluons (jets) [3,4]. 

In the following, we will comment only briefly on jets, in order to point out the 
relation between the two probes, and then concentrate on quarkonium production. 
Quarkonium and jet production are rather well understood in hadron-hadron col- 
lisions, where they are accounted for in terms of perturbative QCD and hadronic 
parton distribution functions [5,6]. In both cases, the initially formed state (QQ : q 
or g) is in general coloured, and it has an intrinsic mass or momentum scale much 
larger than Aqcd- For jets, this is also the state to be used as probe, since the 
behaviour of a fast colour charge passing through confined matter differs from that 
in a deconfined medium [7] . In confined matter, the colour charge loses energy as it 
passes from one hadron to the next through the "interhadronic" vacuum, and the 
energy loss is determined by the string tension a acting on the colour charge as it 
leaves the field of a hadron [8] . In a hot deconfined medium, the crucial quantities 
are the colour screening radius and the mean free path; these determine with how 
many other charges the passing colour charge interacts and how much energy it 
loses per unit length [9,10]. The fate of a colour charge in the transition region 
between these two limits is still quite uncertain. 

For fast quarkonia, the situation is similar; they will pass through the medium 
while still in a coloured state [8] , and hence they can be used as probe in the same 
way as jets. In addition, however, we can consider slow quarkonia, which have be- 



come full physical resonances within a hadronic volume around the QQ formation 
point and thus traverse the medium as colour singlets. Since the intrinsic spatial 
scales of J/ip and T, determined by the heavy quark masses and the binding ener- 
gies, nevertheless remain much smaller than the hadronic size Aq^ d , they interact 
only with the partons within a big, light hadron and not with the hadron as a whole. 
They are thus able to probe the partonic state of any medium. In particular, they 
are essentially unaffected by the soft gluons which make up confined matter, while 
the hard gluons present in a QGP will dissociate them [2]. 

For both quarkonia and jets, thermal production in the expected temperature 
range (T < 0.5 GeV) is excluded by the mass or momentum scales involved; we can 
therefore be quite sure that such signals were produced prior to QGP formation. 
They will also not reach an equilibrium with later stages of the medium. Hard jets 
and fast quarkonia require too much of an energy loss for this, while slow quarkonia, 
as noted, are either dissociated or not affected by the medium. 

For both proposed probes, initial state nuclear effects can occur before QGP 
formation. Primary quarks and gluons may undergo multiple scattering or expe- 
rience shadowing in the nucleus before they interact to form a QQ pair or a hard 
transverse parton. These effects have to be understood and taken into account 
before any QGP analysis [11,12]. It is therefore necessary to study them in pro- 
cesses which are not effected by the subsequent medium, such as the production of 
hard direct photons [13] or of high mass Drell-Yan dileptons [14]. In these cases, 
we have only annihilation or bremsstrahlung of the incident partons; the resulting 
electromagnetic signal leaves the system unaffected by any subsequent medium and 
its evolution. If such reactions show nuclear effects, they are presumably due to 
initial state phenomena. 

After these more general remarks, we will now consider in detail the use of 
quarkonium production as a probe for deconfinement in dense strongly interacting 
matter. We concentrate on quarkonium for several reasons. J/ip suppression was 
predicted [1] to be the consequence of QGP formation, and such a suppression was 
subsequently indeed observed in high energy nuclear collisions at the CERN-SPS 
[15]. This triggered an intensive study of possible alternative origins of such a 
suppression. Hence the analysis necessary to establish an unambiguous probe for 
deconfinement has been carried much further here than for jets and can provide 
a good illustration of what needs to be done before drawing any conclusions. In 
particular, as noted at the beginning, we must understand theoretically and exper- 
imentally the dynamics of the process to be used as probe and 

- how it is influenced by initial state nuclear effects, 

- how it reacts to confined matter, 

- how it reacts to deconfined matter, and 

- how it reacts to non-equilibrium systems. 



In section 2, we will therefore outline the theory of quarkonium production and 
compare it to present data [6]. Section 3 will deal with quarkonium production 
in hadron-nucleus collisions. In particular, we will see how the production of fast 
quarkonia provides us with information on the energy loss of a colour charge in 
confined matter and on gluon shadowing. Next, in section 4, we will discuss the 
conceptual basis of quarkonia as deconfinement probe and present the results of the 
heavy quark theory for quarkonium-hadron interactions. We then show how to test 
experimentally the resulting transparency of confined matter to slow J/ifj's. Section 
5 will bring a comprehensive comparison of the effect of confined and deconfined 
media on quarkonium production; here we will recover the melting of the J/ift in 
a QGP [1,16] on a microscopic level. In Section 6, we remove initial state nuclear 
modifications and study the effect of pre-equilibrium deconfinement. Finally, in 
section 7, we give a a brief summary and an assessment of what we have learned 
from present data. 

2. Quarkonium Production in Hadron-Hadron Collisions 

In this section, we shall sketch the basic dynamics of quarkonium production 
in hadron-hadron collisions; for a summary, see e.g. [6]. We shall speak about 
charmonium states; but everything said also holds for bottonium. 

2.1 Colour Evaporation 

The first stage of charmonium formation is the production of a cc pair; because 
of the large quark mass, this process can be described by perturbative QCD (Fig. 
2.1). A parton from the projectile interacts with one from the target; the parton 
distributions within the hadrons are determined e.g. by deep inelastic lepton-hadron 
scattering. Initially, the cc pair will in general be in a colour octet state. It subse- 
quently neutralises its colour and binds to a physical resonance, such as J/ip, Xc or 
ifi'. Colour neutralisation occurs by interaction with the surrounding colour field; 
this and the subsequent resonance binding are presumably of non-perturbative na- 
ture ("colour evaporation" [17]). In the evaporation process, the cc pair can either 
combine with light quarks to form open charm mesons (D and D) or bind with 
each other to form a charmonium state. 

The basic quantity in this description is the total sub-threshold charm cross 
section, obtained by integrating the perturbative cc production over the mass inter- 
val from 2m c to 2m d- At high energy, the dominant part of a c5 comes from gluon 
fusion (Fig. 2.1a), so that we have 




(2.1) 



with g p (x) and gt(x) denoting the gluon densities in projectile and target, respec- 
tively, and a the gg — > cc cross section. In pion-nucleon collisions, there are also 
significant quark-antiquark contributions (Fig. 2.1b), which become dominant at 
low energies. The essential prediction of the colour evaporation model is that the 
production cross section of any charmonium state i is given by 

0i(s) = ft &cc(s), (2.2) 

where fi is a constant which for the time being has to be determined empirically. 
In other words, the energy dependence of any charmonium production cross section 
is predicted to be that of the perturbatively calculated sub-threshold charm cross 
section. As a consequence, the production ratios of different charmonium states 

44 = T = const - ( 2 - 3 ) 

are predicted to be energy-independent. - We note that in the generalised colour 
evaporation model [6], only a part of the total subthreshold cross section a c5 goes 
into charmonium formation. In accord with perturbative open charm calcualtions, 
the remainder (more than 50 %) leads to DD production, with the missing energy 
obtained by interaction with the colour field. 



Fig. 2.1: Production of a cc pair by gluon fusion (a) and qq annihilation (b) 

2.2 Quarkonium Production: Theory and Data 

The predictions of the colour evaporation model have recently been compared 
in a comprehensive survey [6] to the available data, using parton distribution func- 
tions [18,19] which take into account the new HERA results. In Figs. 2.2 and 2.3, 
we see that the energy-dependence is well described for both J/ip and T production; 
for J/ip production, the normalisation coefficient is /jy^=0.025. The T results are 



obtained for the sum of T, T' and T" decaying into dimuons, with Bfc = 1.6 x 10 -3 
for the normalisation; here the branching ratios cannot be directly removed. In the 
fixed target/ISR energy range, the results from the two different sets of parton dis- 
tributions coincide; for the J/ip at LHC energies, there is some spread due to scale 
uncertainties in the parton distributions, which hopefully can be removed by more 
precise DIS data. For the T production, we have already now data up to 1.8 TeV, 
and in Fig. 2.3 they are seen to agree very well with the prediction obtained using 
the "low energy" value Bfx = 1.6 x 10 -3 . Previous phenomenological fits [20] had 
led to much smaller rates; they are included (CR) in Figs. 2.2 and 2.3. 

In Figs. 2.4 and 2.5, the predicted energy-independence of production ratios 
is found to hold as well, again up to Tevatron energy. Here it should be noted 
that the CDF data for the ratio ip'/(J/ip) are taken at large transverse momenta 
(5 < pr < 15 GeV), while the lower energy data are integrated over px, with the 
low pr region dominant. Hence colour evaporation appears to proceed in the same 
way at both small and large pr- 

The colour evaporation model does not determine the relative production rates 
of the different states. In order to obtain them, the colour evaporation process has 
to be specified in more detail. As an example, we consider the ratios of the different 
I = states shown in Figs. 2.4 and 2.5. Assume that the initial colour octet state 
first neutralises its colour by interaction with the surrounding colour field, producing 
a colour singlet cc state. The relative weights for J/ip and ip' production can then 
be expressed [21] in terms of the corresponding masses and the squared charmonium 
wave functions at the origin, 

Here ip denotes the directly produced IS cc state, in contrast to the experimentally 
observed J/ip, 40% of which originates from radiative Xc decays [22,23]. The wave 
functions at the origin can in turn be related to the dilepton decay widths T ee ~ 
(i? 2 (0)/M 2 ) [24], giving 



Inserting the measured values for masses and decay widths, we find 

— y-rt- ~ 0.24 . 2.6 

To compare this to the measured value of {a{^)') / a{^))), we have to remove the Xc 
contributions from the experimental ratio, 

«_m = r i l r «_m i (2 7) 



With the experimental values a(tp')/a(J/ip) ~ 0.14 (see Figs. 2.2 and 2.3) and 
(a x /(jj) ~ 0.4 [22,23], this yields a(ip')/a(ip) ~ 0.23, in good agreement with the 
theoretical result (2.6). We thus find that the projection of the colour singlet cc 
state onto J/ip and ip' correctly describes their production ratios at all energies and 
transverse momenta. 

The predictions for direct bottonium production ratios corresponding to Eq. 
(2.6) are 

- 0.36 ; ^1 * 0.27 . (2.8) 
a(T) (t(T) v ' 

Since the contributions from indirect production through radiative Xb decay are 
not yet known and there is also feeddown from higher S-states, a quantitative 
comparison is not possible here. Nevertheless, the predicted values differ by less 
than 50 % from the data and hence appear reasonable. 

3. Quarkonium Evolution and Hadron-Nucleus Collisions 

Quarkonium production in a (finite) medium depends crucially on the evo- 
lution stage in which the QQ pair is during its passage. Here we determine the 
different stages attainable in hadron-nucleus collisions and discuss the effects which 
the nuclear environment has on them. 

3.1 Parton Fusion and Colour Neutralisation in Nuclear Matter 

The first stage of charmonium production by hadronic interactions is, as we 
saw, parton fusion resulting in a generally coloured cc pair. Subsequently, the cc 
becomes colour neutral by emission or absorption of gluons in the colour field of 
the interaction, and eventually this small colour singlet state "expands" to form a 
physical charmonium resonance. If all or part of this evolution takes place inside a 
nucleus, a number of effects can modify the production: 

- the effective distribution of the partons which fuse to form the cc can be mod- 
ified by the nuclear environment (EMC effect); 

- the colour octet cc state will interact with the nuclear medium strongly and 
without knowledge of its final charmonium state (nuclear shadowing, cc energy 
loss); 

- the physical charmonium states will interact with the nuclear medium accord- 
ing to their different cross sections on hadrons (absorption). 

Since the evolution is temporal, with finite time scales, fast and slow charmonia 
(in the nuclear rest frame) will have quite different fates. In this section, we shall 
determine the kinematic regimes for the different evolution stages [8,25]; in the 
subsequent sections, we will then discuss the dominant effects in each stage. 

To estimate how long the colour octet cc state will live, we first note that the 
J/ip wave function keeps the charmed quarks close to their mass shells; we thus have 



p2 ^ m 2 ^ gee pjg 3 ^ £ or notation). The intermediate quark with four-momentum 
k = p + q is off-shell by an amount 

A = (p + q) 2 -m 2 c = 2pq, (3.1) 

where q is the momentum of the third (colour neutralising) gluon; in the spirit of the 
parton model, we keep all gluons on-shell. In the low pr region, in which we are here 
primarily interested, this third gluon can be arbitrarily soft; colour neutralisation 
could even involve several gluons. In any case, the colour neutralisation process 
cannot really be calculated in a purely perturbative framework. We shall here 
nevertheless keep the structure of Fig. 3.1, taking it to be a phenomenological 
extension into the non-perturbative soft-gluon regime. A justification of such a 
procedure might come from recently developed resummation methods [26]. 



Fig. 3.1: J/ifj production through gluon fusion 

The proper life-time of the virtual coloured state is now given by the uncer- 
tainty principle as 

rg^lAI" 1 / 2 . (3.2) 
In the rest frame of the J ftp, the life-time of the colour octet, 

r 8 ~ (2m c?0 )- 1/2 , (3.3) 

is determined by the energy q of the third gluon. It is clear from Eqs. (3.1) and 
(3.2) that the colour octet state can propagate over long distances only if the third 
gluon is soft enough. In a confined medium, gluons cannot propagate over distances 
larger than about 1 fm; hence the low energy cut-off is go — Aqcd — 0.2 GeV. This 
leads to 

r 8 ~ 0.25 fm (3.4) 



for the colour neutralisation time. In the rest frame of the nuclear target A, the cc 
travels in this time a distance 



where Pa is the (lab) momentum of the cc and M its mass. From Eq. (3.5) it is 
clear that in hadron-nucleus collisions sufficiently fast cc pairs will still be coloured 
when they leave the nuclear medium. 

The average path length for a cc produced in a h — A collision is (3/4)Ra — 
0.86A 1 / 3 fm; for heavy nuclei, such as Pb (A = 208), this becomes about 5 fm. 
Hence charmonium states of lab momenta Pa > 60 GeV have passed the nucleus 
as colour octets. To relate this to the kinematic variables generally used in h — A 
experiments, we transform the lab momentum Pa to the momentum P measured in 
the center of mass of a hadron-nucleon collision. The Feynman variable xp is then 
defined as Xp = P/P-max, with P m ax{s) denoting the maximum cms momentum 
possible for a charmonium state produced in a hadron-nucleon collision of cms 
energy y/s. From this, we obtain [8,25] that for y^s > 20 GeV and xp > 0, the 
measured charmonium states have traversed most of the nucleus as colour octets. 
This kinematic range covers essentially all high statistics charmonium production 
data taken in h — A collisions [27,28]. These experiments therefore study the passage 
of a colour octet cc pair through nuclear matter; as such, they do not provide any 
information on the interaction of fully formed physical J/ip or ip' with nucleons. 

An immediate consequence of this is that any nuclear effects observed for char- 
monium production in h — A collisions in the quoted kinematic range must be 
the same for J/ip and for ip' production. This conclusion is indeed well satisfied 
experimentally (Fig. 3.2). 

We have so far considered the colour regime, i.e., the range in which the cc 
pair traverses the entire nucleus as a colour octet. When the distance do is less 
than 1-1.5 fm, the cc is colourless when it passes from the nucleon on which it 
was produced to the "next" nucleon inside the nucleus. At yfs — 20 GeV, this 
occurs for xp < — 0.2. At this point, the cc is still a generic small colour singlet 
state and not yet a specific physical resonance. The other extreme to the colour 
regime is the range in which the cc has become a fully formed physical resonance 
{J/ip, ip') before it leaves the environment of the nucleon at which it was formed. 
Since the different states have different sizes and hence different formation times, 
the resonance regime will be state-specific, i.e., at a given collision energy, it will 
cover a larger ^-region for the J/ip than for the ip' '. The distance d r (i) travelled 
by the cc before becoming a physical resonance state i is given by 




(3.5) 




(3.6) 



where r r (i) denotes the resonance formation time. Potential theory [29] leads to 
the estimates t r (J/ij)) ~ 0.35 fm and t r {^') ~ 1 fm. For a collision energy y/s ~ 20 
GeV, this means that a J/ip has become a fully formed resonance for xf < — 0.5. 
The ip' is still nascent even for xf = — 1; although it is a colour singlet on its 
passage through the nucleus if xf < — 0.2, it has not yet reached its full physical 
size when it leaves the nucleus even for xp = — 1. 

In the transition regime between these two extremes, the cc experiences colour 
interactions on part of its passage, while it traverses the medium as a small colour 
singlet on the remaining part. In Table 1, we have summarised the different kine- 
matic ranges for proton beam incident on a Pb target at a cms nucleon-nucleon 
collision energy y/s = 20 GeV, for (IS) J/ip and (2S) ip' production. In the follow- 
ing sections, we survey the possible nuclear effects which can arise in the different 
kinematic regimes. 



Regime 






Colour Regime 


X F > 


X F > 


Transition Regime 


-0.5 < x F < 


-1 < X F < 


Resonance Regime 


-1 < x F < -0.5 





Table 1: 

Charmonium Production in p — A Collisions at ^/s = 20 GeV 



3.2 Quantum-Mechanical Interference and Nuclear Shadowing 

We now address nuclear modifications of charmonium production due to what 
at first sight appears to be a change in the effective parton distribution function 
gt{xt) in nuclei, compared to that in nucleons. The fraction of the target momentum 
carried by the corresponding parton in the elementary interactions shown in Fig. 
2.1 is given by 

x t = ^(4Mys) + x 2 F ±x F ), (3.7) 

where M is the mass of the cc system; the ± signs hold for xf > and xf < 0, 
respectively. For cc production at xp > and y/s > 20, x t < 0.15. In this x t 



region, deep inelastic scattering experiments on nuclei find a lower quark density 
than on nucleons (nuclear shadowing [30]), as illustrated in Fig. 3.3 (the solid line 
is phenomenological fit [31]. Such nuclear shadowing can, however, not be inter- 
preted as an intrinsic change of the parton distribution in nuclei. As such, it would 
be applicable in factorisable form to all hard processes; but for Drell-Yan dilepton 
production by p — A collisions, which involves nuclear quark distribution in the 
same Xt region, little or no modification is observed [28]. The appearent puzzle is 
resolved by noting that nuclear shadowing is due to quantum-mechanical interfer- 
ence, similar to the well-known Landau- Pomeranchuk effect [32] . This interference 
depends on the specific process and hence leads to different nuclear modifications for 
charmonium hadro- and photoproduction, for Drell-Yan dilepton production and 
for deep inelastic lepton-hadron scattering [12]. For a first discussion of quantum 
interference effects in deep inelastic scattering, see [33]. 

We had seen above in Eq. (3.5) that sufficiently fast cc pairs traverse the entire 
nucleus in a virtual coloured state. In this case, the interactions of the pair with 
different nucleons cannot be added incoherently and factorisation breaks down even 
for dynamically uncorrelated nucleons: quantum mechanical interference can now 
lead to nuclear modifications even though both the parton distribution and the 
elementary parton interaction amplitude are unchanged. Such interference effects 
set in when the coherence length d c = (1/mxt), over which the cc is in an off-shell 
virtual state, becomes larger than the internucleonic distance do = n ~ 1.8 fm, 
with no = 0.17 fm -3 denoting standard nuclear density and m the nucleon mass. 
This distance has to be compared to the mean free path A§ of the virtual colour 
state in nuclear matter, 

A 8 = (a 8 n )-\ (3.8) 

where ag is the cross section for the interaction of the colour octet cc with nucleons. 
To estimate the size of cr 8 , we assume that the overall colour neutrality required for 
the propagation of the coloured cc through the nucleus is provided by a comoving 
light qq pair. The interaction cross section of the QQqq system is then defined 
by the light qq pair, so that its size can be as large as that of light qq mesons 
(20 - 30 mb). This is much larger than that of a physical J/ip; it is now not 
the size of the cc, but its colour which determines the interaction. The mean free 
path corresponding to such a cross section becomes Ag ^ 2 — 3 fm, which is about 

— 1/3 

the size of the internucleonic distance do = n ~ 1.8 fm, so that charmonium 
production in p — A collisions for xp > and yfs > 20 GeV satisfies the "shadowing 
relation" <i — A 8 <g; d c [12]. In this regime, the interactions of the virtual cc 
with successive nucleons interfere destructively with each other; it is the regime 
which for electromagnetic interactions results in the Landau-Pomeranchuk effect. 
The production cross section a{pA — > cc) on nuclei is here therefore less than the 



incoherent result Aa(pp — > cc). 

The basic quantity studied experimentally is the ratio 

a(pA -> cc) 

Ra ' p ~ Aa{pp - cc) ' (3 ' 9) 

it would be unity in the absence of any nuclear effects. It is observed, however, 
that Ra/ p decreases with increasing cc momentum in the nuclear rest frame (i.e., 
whenever d c increases), as well as with increasing A. Rather than in terms of 
d c , Ra/p is generally studied in terms of the fractional target parton momentum 
Xt = (m<i c ) _1 . Decreasing xt thus increases the path length for the colour octet 
inside the nucleus, and a longer path of the virtual state leads to more destructive 
interference. When d c reaches the nuclear diameter 2Ra, we have the maximum 
possible interference, so that the suppression saturates for x t < xf at ~ (2t7iRa)~ 1 , 
with the approximate limit [12] 

<yo * ^- (3.io) 

A more realistic form including geometric effects is given by 

R% p {xr)^l-e- x ^ R \ (3.11) 

which reduces to Eq. (3.10) for X/Ra 1. The available data fall into the region 
< Ra/^s = noagRA ^ 0.7, so that the measured nuclear reduction lies between 
0.5 and 1. In this region, the relation e~ x ~ 1 — e 1 ^ is quite well satisfied, so that 

R A/ P { X T) - e~ noasRA (3.12) 

also should also describe the data. The form (3.12) was in fact obtained as an 
empirical fit [34]; we see here, however, that it is the result of quantum-mechanical 
interference effects of the virtual cc and is not related to the interaction of a physical 
J/ip with nucleons, as assumed there. 

For d c < do, the scattering on different nucleons becomes incoherent, giving 
R A h / p = 1) with x t > x™ c ~ (m<io) _1 - For A = 200, the two limits are x s t at = 0.015 
and xf 10 = 0.12. Qualitatively, this is the type of behaviour found in charmonium 
production for xp not too large [11]. The specific functional form of the suppression 
between the two extremes can only be estimated phenomenologically. Experimental 
results have often been parametrised as Ra/ p — A a_1 ; our interference consider- 
ations imply that the coefficient a depends on xt- An analysis of charmonium 
production in p — A collisions [11] suggests a simple linear fit in In A and lnx t: 



Ra/p-^-c Hxtmd )lnA}. 



(3.13) 



The scale mdo ~ 0.12 is chosen to make R 8 a/ p = 1 in the incoherence regime; the 
constant c has to be chosen such as to give the correction saturation value. For a 
fit of p — A data to a similar form, see [8] . 

3.3 Momentum Loss of a Colour Charge in Nuclear Matter 

The suppression R s J^/ p cannot be directly compared to data, however, since 
shadowing is not the only effect suffered by the coloured cc state on its path through 
the nucleus. The cc is formed in the collision of the projectile with one of the target 
nucleons. To leave the nuclear medium, it has to traverse the remaining part of the 
nucleus, in which it can interact with other nucleons. Getting from one nucleon 
to the next means passing regions of internuclear physical vacuum, which does not 
support colour charges. To estimate the effect of this passage, we imagine the cc 
to stretch a string from its formation region to the next nearest nucleon, then from 
this to the one after that, and so on. If the pair intially had a cms momentum P, 
then this will be shifted to 



P' = P-K 8 (L A -d ), (3.14) 

where Kg — (9/4) k ~ (9/4) GeV/fm is the string tension of the colour octet, k ~ 
1 GeV/fm that of the fundamental triplet [8]; La is the total path length of the 
colour state. As consequence of this momentum loss, a cc pair observed at a given 
x f must have been originally produced at a higher value xf/5, with 

P- KB (L A -do)\ (3J5) 



The resulting normalised ^-distribution Fa(xf) = (l/A)(da h A^j/^/dxF) for 
J/ip production in h — A collisions thus becomes 

Fa(xf) = WoF (xf) + (1 - Wo) ( F ° ( ^ F/a) ) 9(1 - x F /S). (3.16) 

Here F (xf) is the ^-distribution in h — p collisions and 

Wo ~ exp{-(L A - z )/d 8 } (3.17) 

gives the probability that the cc pair will not undergo any scattering on its way 
out of the target; thus Wo = 1 implies no nuclear J/ip suppression. The second 
term in Eq. (3.16) gives the effect of the n coherent scatterings in the medium, 
with the resulting shift in xp- The factor 1/8 assures that the distribution remains 
normalised. Note that the momentum loss of cc pair, for xf > as described here, 
does not imply any integrated J/ip suppresion: the production is simply shifted from 



larger to smaller xp. Dividing Eq. (3.16) by the ^-distribution for h-p collisions, 
Fq(xf), gives us the ratio of h — A to h — p distributions, 



For a nucleus of A = 200 and the other parameters as above, we get the dependence 
on xf shown in Fig. 3.4 for three different beam energies. We note that on nuclei, 
compared to nucleons as target, J/ip production is shifted to lower xf, essentially 
because of the momentum degradation of the colour octet in passing through the 
medium. This momentum loss saturates at AP max ~ k 8 (Ra — zo) — 10 GeV and is 
thus bounded. As a consequence, the relative shift in xf decreases with increasing 
beam energy, and the effect of the medium disappears for y/s — > oo, when R™j p — 1 
for < xp < 1. 



Fig. 3.4: The effect of colour octet interactions on charmonium production in p — A 
collisions at three different energies 

In addition to the non-perturbative nuclear modification of fast quarkonium 
production in nuclei, the cc pair can undergo hard interactions with the quarks and 
gluons within each of the nucleons it passes through. We neglect this here, since the 
density of sufficiently hard scatterings in nuclei is too low to compete with string 
stretching, as will become evident in Section 3.4. This is changed, however, in a 
quark-gluon plasma, and the difference in parton hardness for confined and decon- 
fined media is in fact crucial for the use of quarkonia as confinement/deconfinement 
probe. 

The observed suppression of quarkonium production in h — A collisions will be 
due to both the effects just discussed, nuclear shadowing and momentum loss of 




(3.18). 



the coloured cc pair. We thus have to compare experimental results to the product 

R A/p (x F ) = R s A ) p (x F ) x R% p {x F ) (3.19) 

of the two mechanisms, as given by Eqs. (3.12) and (3.18). The energy loss is 
effective mainly at larger x F ; since the integrated cross section, on the other hand, 
is determined by the region around x F ~ 0, its suppression is dominantly due to 
nuclear shadowing. Hence Eq. (3.11) should determine the main A-dependence of 
charmonium suppression in p — A collisions at Xf > 0, and as a comparison to the 
equivalent form (3.12) shows [34], it indeed does. 

This form, but with a slightly different parametrisation of shadowing [11], has 
been compared to the available data [27,28] in [8]. In spite of the simplistic fit to 
shadowing and some kinematic approximations at small x F , the result (Fig. 3.5) is 
seen to reproduce the observed suppression quite well. 

Quarkonium production in hadron-nucleus collisions, as experimentally stud- 
ied up to now, appears thus as theoretically quite well understood, based on the 
interaction of a fast colour octet cc pair with nuclear matter. The dominant mech- 
anisms for this interaction are the destructive quantum-mechanical interference of 
the scattering amplitudes of the cc on different nucleons, and the momentum loss of 
the cc as it traverses the physical vacuum between successive nucleon interactions. 

3.4 Charmonium Interactions in Nuclear Matter 

If we want to know the effect of nuclear matter on a fully formed physical 
J/ip, then - as shown in Table 1 - we have to study the production of J/^'s slow 
in the nuclear rest frame. The study of such slow charmonia has up to now been 
essentially impossible. It requires the detection of slow decay dileptons, and for 
this the abundance of slow hadrons constitutes an overwhelming background. In 
the case of fast dileptons, a hadron absorber can eliminate these, and hence all 
p — A studies were so far restricted to dilepton pairs of more than 20 GeV in the 
rest frame of the nuclear target. The advent of the P6-beam at the CERN-SPS has 
removed this constraint. With the P6-beam incident on a hydrogen (or deuterium) 
target, the nuclear rest frame moves with a lab momentum of 160 GeV. Hence 
now those charmonia (and their decay dileptons) which are slow in the nuclear rest 
frame are very fast in the lab system and will thus pass the hadron absorber; such 
experiments can therefore provide the cross section for the break-up interaction of 
physical J/ip : s on nucleons. To avoid confusion, we shall continue to define x F also 
in such Pb — h reactions as positive for cms momenta in the direction of the hadron; 
with this terminology, Pb — h collisions provide us with information about the so 
far unknown region of negative x F . Estimates [8,25] show that for yt~s > 17 GeV 
and x F < —0.2, a cc is colourless before it reaches the next nucleon, so that colour 
interactions no longer enter. 



We now parametrise the survival probability for state i (i = J/ip, Xo VO as 



Si = exp{—n (JiL} , (3.20) 

where no is as before normal nuclear density, L the average path length of the 
charmomium state i in the nucleus, and Oi its absorption cross section of the state 
% in nuclear matter. If the charmonium state is fully formed before it leaves the 
range of the nucleon at which it was produced, <7j is simply the charmonium-nucleon 
cross section. If it is not yet fully formed, the effective cross section will be smaller, 
vanishing in the colour transparency limit of a pointlike colour singlet [35]. For 
illustration, we can therefore use a simplistic parametrisation of the cross section 
as function of the distance d the state has travelled before leaving the nucleus 
[36,37], 

a i {d) = a i (i-\ , (3.21) 



where Oi is the fully developed cross section and L = [(Pa/M^t^ — 1] the effective 
distance travelled until full resonance formation. Eq. (3.21) holds for d < L; for 
d>L, o~i(d) = Oi. Using this parametrisation, Eq. (3.20) is replaced by 

2 - 

Si = exp{-n cxi[L - g£;]}; (3.22) 

this reduces to Eq. (3.20) for momenta high enough to make L = 0. To determine 
the actual survival probabilities, we now need the values Oi for the fully formed 
resonances colliding with nucleons. That is a topic in its own right and will be 
taken up in the next chapter. Before turning to this question, we want to connect 
the survival probability (3.20) to the measured suppression ratio Ra/ p - 

3.5 EMC Effect Modifications 

The survival probabilities Si are related to the p — A and p — p induced char- 
monium production cross sections through 

(dof A /dx F ) 9a{x) 

A( d *r/dx F ) = ^m Si(xf) = R ^ S ^ ( 3 - 23 ) 

where qa(x) and qa(x) are the parton distribution functions in nuclear and proton 
target, respectively. The fractional parton momentum x is now, with xp < 0, 
expressed by 

x= l -{^{xl + {AMf/s))-x F ) (3.24) 

in terms of the variables xf and s. At the energy possible for P6-beam experiments 
(y/s = 17.4 GeV), and for xp < 0, we have x > 0.15. We are thus above the region 



in which quantum-mechanical coherence effects (nuclear shadowing) can play a role; 
the reaction now is indeed an incoherent sum of interactions with different nucleons. 
In deep inelastic scattering on nuclear targets there is for x 2 > 0.15 a modification 
in comparison to the same process on nucleons (EMC effect [38]; it can now be 
interpreted as a genuine change of the parton distribution function and applied to 
other hard processes. The resulting pattern for quark distributions q(x), 



R A/P M-—^ ( 3 - 25 ) 



is illustrated in Fig. 3.6 for M = 3.1 GeV and ^fs = 17 A GeV, making use of Eq. 
(3.24) to convert the Bjorken variable x to xp. 

Although the EMC effect has so far been observed only for quarks, it is to be 
expected that gluons will exhibit a similar behaviour. This is one reason why the 
initial state factor R^/ p (x) will introduce an further xp variation in addition to that 
coming from Si(xp). A second reason was already indirectly mentioned above: with 
increasing \xf\-, charmonium production is more and more due to quark-antiquark 
annihilation rather than to gluon fusion. The two contributions become approxi- 
mately equal around \xp\ = 0.5, and as \xp\ — > 1, the qq contribution is dominant. 
We shall here assume that the gluon distributions behave similarly and use the 
quark form of Ra/ p m the whole region —1.0 < xp < 0. Eqs. (3.22)/(3.23) and the 
Rjifp °f Fig- (3.6) then predict the behaviour of the measured cross section ratio 
in the region of negative xp. 

4. The Theory of Quarkonium-Hadron Interactions 

In this section, the cross sections for the inelastic interaction between quarkonia 
and light hadrons will be calculated in short distance QCD, based on the small 
radii and the large binding energies of the lowest QQ resonances. A break-up of 
such states requires hard gluons, but the gluons confined to hadrons slow in the 
quarkonium rest frame are in general very soft, so that the resulting cross sections 
remain very small until quite high collision energies. 

4.1 The Short-Distance Analysis of Quarkonium-Hadron Interactions 

The interaction of quarkonia with ordinary light hadrons plays an important 
role both in the dynamics and in the thermodynamics of strong interaction physics. 
For QCD dynamics, it is important since the small quarkonium size probes the 
short-distance aspects of the big light hadrons and thus makes a parton-based 
calculation of the overall cross section possible [39-43]. In QCD thermodynamics, 
quarkonia can be used as a probe for deconfinement [1,2], provided their interaction 
in dense confined matter can be distinguished from that in a quark-gluon plasma. 



We begin this chapter with a brief summary of the QCD analysis of quarko- 
nium interactions with light hadrons. It concludes that the small size of quarkonia, 
combined with the rather large mass gap to open charm or beauty, strongly in- 
hibits their break-up by low energy collisions with light hadrons [40]; the total 
quarkonium-hadron cross section attains a constant asymptotic value only at very 
high energies, compared to corresponding cross sections for the interaction between 
light hadrons. This slowly rising form of the cross section is derived from an oper- 
ator product expansion with ensuing sum rules and becomes quite transparent in 
parton language. 

The QCD analysis of quarkonium interactions applies to heavy and strongly 
bound quark-antiquark states [40] ; therefore we here restrict ourselves to the lowest 
cc and bb vector states J/ip and T, which we denote generically by <E>, following the 
notation of [40]. For such states, both the masses tjiq of the constituent quarks 
and the binding energies eo(<&) — (2M(g g ) — M$) are much larger than the typical 
scale Aqcd for non-perturbative interactions; here (Qq) denotes the lowest open 
charm or beauty state. In $ — h interactions, as well as in ^-photoproduction, 
7/1 — > we thus only probe a small spatial region of the light hadron h; these 
processes are much like deep-inelastic lepton-hadron scattering, with large tuq and 
e in place of the large virtual photon mass \f—q 2 . As a result, the calculation of 
photoproduction and of absorptive $ — h interactions can be carried out in the short- 
distance formalism of QCD. Just like deep-inelastic leptoproduction, these reactions 
probe the parton structure of the light hadron, and so the corresponding cross 
sections can be calculated in terms of parton interactions and structure functions. 

In the following, we shall first sketch the theoretical basis which allows quarko- 
nium interactions with light hadrons to be treated by the same techniques as used 
in deep-inelastic lepton-hadron scattering or in the photoproduction of charm. We 
show the derivation for the sum rules which relate the absorptive $ — h cross sec- 
tion to hadronic gluon structure functions [39,40]. This relation given, we calculate 
explicitly the energy dependence of the cross section. Readers only interested in 
this behaviour can therefore go immediately to Eq. (4.24). 

Consider the amplitude for forward scattering of a virtual photon on a nucleon, 



In the now standard application of QCD to deep-inelastic scattering one exploits the 
fact that at large spacelike photon momenta q the amplitude is dominated by small 
distances of order 1/^—q 2 (Fig. 4.1a). The Wilson operator product expansion 
then allows the evaluation of the amplitude at the unphysical point pq — > 0, where 
p is the four-momentum of the nucleon. Since the imaginary part of the amplitude 
(4.1) is proportional to the experimentally observed structure functions of deep- 
inelastic scattering, the use of dispersion relations relates the value of the amplitude 




(4.1) 



at pq — > point to the integrals over the structure functions, leading to a set of 
dispersion sum rules [44]. The parton model can be considered then as a particularly 
useful approach satisfying these sum rules. 

In the case J M = Q^^Q, i.e., when vector electromagnetic current in Eq. (4.1) 
is that of a heavy quark-antiquark pair, large momenta q are not needed to justify 
the use of perturbative methods. Even if q ~ 0, the small space-time scale of x 
is set by the mass of the charmed quark, and the characteristic distances which 
are important in the correlator (4.1) are of the order of l/2mg (Fig. 4.1b). In 
[41,42], this observation was used to derive sum rules for charm photoproduction 
in a manner quite similar to that used for deep-inelastic scattering. 

In the interaction of quarkonium with light hadrons, again the small space scale 
is set by the mass of the heavy quark, and the characteristic distances involved are of 
the order of quarkonium size, i.e., smaller than the non-perturbative hadronic scale 
Aq£ D . Moreover, since heavy quarkonium and light hadrons do not have quarks in 
common, the only allowed exchanges are purely gluonic. However, the smallness of 
spatial size is not enough to justify the use of perturbative expansion [40]. Unlike in 
the case of $-photoproduction, heavy quark lines now appear in the initial and final 
states (see Fig. 4.1c), so that the QQ state can emit and absorb gluons at points 
along its world line widely separated in time. These gluons must be hard enough to 
interact with a compact colour singlet state (colour screening leads to a decoupling 
of soft gluons with the wavelengths larger than the size of the $); however, the 
interactions among the gluons can be soft and nonperturbative. We thus have to 
assure that the process is compact also in time. Since the absorption or emission of 
a gluon turns a colour singlet quarkonium state into a colour octet, the scale which 
regularizes the time correlation of such processes is by the quantum-mechanical 
uncertainty principle just the mass difference between the colour-octet and colour- 
singlet states of quarkonium: r c ~ l/(es — ei). The perturbative Coulomb-like piece 
of the heavy quark-antiquark interaction 

V k (r) = -g>£- (4-2) 

is attractive in the colour singlet (k = 1) and repulsive in the colour-octet (k = 8) 
state; in SU(N) gauge theory 

q 2 N 2 -l 

* = -h— < 4 ' 3 «> 

To leading order in 1/N, the mass gap between the singlet and octet states is 
therefore just the binding energy of the heavy quarkonium eo, and the characteristic 



correlation time for gluon absorption and emission is 



t c ~ 1/eo- (4.4) 

Although the charm quark is not heavy enough to ensure a pure Coulomb regime 
even for the lowest cc bound states (rj c and J/ip), the mass gap determined from 
the observed value of open charm threshold clearly shows that r c < A-qcd- For 
the T, the interaction is in fact essentially Coulomb-like and the mass gap to open 
beauty is even larger than for charm. One therefore expects to be able to treat 
quarkonium interactions with light hadrons by the same QCD methods that are 
used in deep-inelastic scattering and charm photoproduction. 

We thus use the operator product expansion to compute the amplitude of heavy 
quarkonium interaction with light hadrons, 

F* h = i j d 4 xe i(lx {h\T{J(x)J(0)}\h) = J2cn(Q,m Q ){O n ), (4.5) 

where the set {O n } includes all local gauge invariant operators expressible in terms 
of gluon fields; the matrix elements (O n ) are taken between the initial and final 
light-hadron states. The coefficients c n are computable perturbatively [39] and 
process-independent. As noted above, in deep-inelastic scattering the expansion 
(4.5) is useful only in the vicinity of the point pq — > 0. The same is true for the 
case of quarkonium interaction with light hadrons. As shown in [40], the expansion 
(4.5) can therefore be rewritten as an expansion in the variable 

^ _ PQ_ _ {s-Ml-Mp ^ ( 8 - Ml) 

M$ 2M$ 2M$ 1 ' ' 

where is the mass of the light hadron; the approximate equality becomes valid 
the heavy quark limit. For the lowest IS quarkonium state one then obtains 

oo / A \ n 

F*h = rlel J2 d n(°n)[-) » ( 4 - 7 ) 

where ro and eo are Bohr radius and binding energy of the quarkonium, and the 
sum runs over even values of n to ensure the crossing symmetry of the amplitude. 
The most important coefficients d n were computed in [39] to leading order in g 2 
and 1/N. 

Since the total $ — h cross section a$h is proportional to the imaginary part 
of the amplitude the dispersion integral over A leads to the sum rules 

- / d\ A-V^(A) = r s e 2 d n (O n ) - . (4.8) 



Eq. (4.7) provides only the inelastic intermediate states in the unitarity relation, 
since direct elastic scattering leads to contributions of order r®. Hence the total cross 
section in Eq. (4.8) is due to absorptive interactions only [40], and the integration 
in Eq. (8) starts at a lower limit Ao > Mh- Recalling now the expressions for radius 
and binding energy of IS Coulomb bound states of a heavy quark-antiquark pair, 

ro = ( v ) ^ (4 ' 9) 

(410) 

and using the coefficients d n from [39] , it is possible [40] to rewrite these sum rules 
in the form 

with Ao/eo — 1 in the heavy quark limit. The contents of these sum rules become 
more transparent in terms of the parton model. In parton language, the expectation 
values (O n ) of the operators composed of gluon fields can be expressed as Mellin 
transforms [45] of the gluon structure function of the light hadron, evaluated at the 
scale Q 2 = €q, 

r i 

dx x n - 2 g(x,Q 2 = e 2 Q ). (4.12) 



(On) = r 

Jo 



Defining now 

» = (4-13) 
we can reformulate Eq. (11) to obtain 

/ dy y n - 2 a^ h (X /y) = I(n) [ dx x n ~ 2 g(x,Q 2 = eg), (4.14) 
jo Jo 

with I(n) given by 

Eq. (4.14) relates the $ — h cross section to the gluon structure function. To get a 
first idea of this relation, we neglect the n-dependence of I(n) compared to that of 
(O n ); then we conclude that 



cr^ h (X /x) ~ g(x, Q 2 = el), 



(4.16) 



since all order Mellin transforms of these quantities are equal up to a constant. 
From Eq. (4.16) it is clear that the energy dependence of the $ — h cross section 
is entirely determined by the x— dependence of the gluon structure function. The 
small x behaviour of the structure function governs the high energy form of the cross 
section, and the hard tail of the gluon structure function for x — > 1 determines the 
energy dependence of a$h close to the threshold. 

To obtain relation (4.16), we have neglected the n-dependence of the function 
I(n). Let us now try to find a more accurate solution of the sum rules (4.13). 
We are primarily interested in the energy region not very far from the inelastic 
threshold, i.e., 

(M h + e ) <A< 5GeV, (4.17) 

since we want to calculate in particular the absorption of <fr's in confined hadronic 
matter. In such an environment, the constituents will be hadrons with momenta 
of at most a GeV or two. A usual hadron (jr, p, nucleon) of 5 GeV momentum, 
incident on a J/i/i at rest, leads to y/s ~ 6 GeV, and this corresponds to A ~ 5 
GeV. 

From what we learned above, the energy region corresponding to the range 
(4.17) will be determined by the gluon structure function at values of x not far 
from unity. There the x-dependence of g(x) can be well described by a power law 

g(x)=g 2 (k+1) (l-x) k , (4.18) 

where the function (18) is normalized so that the second moment (4.12) gives the 
fraction g 2 of the light hadron momentum carried by gluons, < 2 >= gi — 0.5. 
This suggests a solution of the type 

<r* h (y) = a(l - y) a , (4.19) 

where a and a are constants to be determined. Substituting (4.18) and (4.19) into 
the sum rule (4.13) and performing the integrations, we find 

r(q + i) = [ ^ 3/2 92 \ fiQ\ 2 fiQn\ r(n + f)r(fc + 2) 

T(n + a) \ m 2 Q )\3J \3g 2 J T(n + 5) T(k + n)' K ' 

We are interested in the region of low to moderate energies; this corresponds to 
relatively large x, to which higher moments are particularly sensitive. Hence for 
the range of n for which Eq. (4.5) is valid [42], n < 8, the essential n-dependence 
is contained in the V- functions. For n > 4, Eq. (4.20) can be solved in closed form 
by using an appropriate approximation for the T-functions. We thus obtain 

a + \ } ~ const. n- fc " 5 / 2 . (4.21) 

r(fc + 2) v ' 



Hence to satisfy the sum rules (4.14), we need 



a = k+ - a = const. -) =f . 4.22 

Therefore the solution of the sum rules (4.13) for moderate energies A takes the 
form 

m n o /2 /16\ 2 /16tt\ 1 r(fc + 2) / A \ fc+5/2 

To be specific, we now consider the J/ ^-nucleon interaction. Setting k = 4 in accord 
with quark counting rules, using g<i ~ 0.5 and expressing the strong coupling g 2 in 
terms of the binding energy eo (Eq. 4.10), we then get from Eq. (4.23) the energy 
dependence of the J/ip N total cross section 

6.5 



a J/i)N (X) ~ 2.5 mbx ( 1 - y J , (4-24) 



with A given by Eq. (6) and Ao — (Mjv + eo)- This cross section rises very slowly 
from threshold, as shown in Fig. 4.2; for Pjv — 5 GeV, it is around 0.1 mb, i.e., 
more than an order of magnitude below its asymptotic value. 

We should note that the high energy cross section of 2.5 mb in Eq. (4.24) is 
calculated in the short-distance formalism of QCD and determined numerically by 
the values of m c and eo- From Eqs. (4.10) and (4.23), it is seen to be proportional 
to l/(mQy/mQeo). For T — N interactions, with m& ~ 4.5 GeV and eo ^ 1.10 GeV, 
we thus have the same form (4.24), but with 

o~yn ^ 0.37 mb (4.25) 

as high energy value. This is somewhat smaller than that obtained from geometric 
arguments [46] and potential theory [47]. 

For sufficiently heavy quarks, the dissociation of quarkonium states by inter- 
action with light hadrons can thus be fully accounted for by short-distance QCD. 
Such perturbative calculations become valid when the space and time scales as- 
sociated with the quarkonium state, tq and £q, are small in comparison to the 
nonperturbative scale Aqq D 

r Q « A qcd> (4.26a) 
tQ « A qcd; (4-266) 

Aqq D is also the characteristic size of the light hadrons. In the heavy quark limit, 
the quarkonium binding becomes Coulombic, and the spatial size tq ~ (a s mQ) _1 



thus is small. The time scale is by the uncertainty relation given as the inverse of 
the binding energy Eq ~ uiq and hence also small. For the charmonium ground 
state J/ip, we have 

rj /4 , ~ 0.2 fm = (1 GeV)" 1 ; E JN> = 2M D - M^> ~ 0.64 GeV. (4.27) 

With Aqcd — 0.2 GeV, the inequalities (4.25) seem already reasonably well satis- 
fied, and also the heavy quark relation Ej^ = (l/m c r^^) is very well fulfilled. We 
therefore expect that the dissociation of J/^'s in hadronic matter will be governed 
by the J/^-hadron break-up cross section as calculated in short-distance QCD. 

Nevertheless, in view of the finite charm quark mass, it makes sense to ask if the 
formalism developed here correctly describes J/ip interactions with light hadrons. 
We will take up this question first theoretically, in the next section, where we 
shall see if non-perturbative interactions can lead to significant contributions to 
J/ip break-up. In Section 4.3, we then see how to verify experimentally the short- 
distance QCD prediction for J/ip break-up by light hadrons. 

4.2 Non-Perturbative Quarkonium Dissociation 

For an isolated J/^-hadron system, non-perturbative interactions can be pic- 
tured most simply as a quark rearrangment [48]. Consider putting a J/ip "into" a 
stationary light hadron; the quarks could then just rearrange their binding pattern 
to give rise to transitions such as J/ip + N^A c + D or J/ip + p — > D + D (Fig. 
4.3). The probability for such a transition can be written as 

-Prearr ~ 

/ d 3 rR(r) |<Mr)| 2 , (4-28) 

where the spatial distribution of the cc bound state is given by the squared wave 
function |^(r)| 2 . The function R(r) in eq. (4.28) describes the resolution capability 
of the colour field inside the light hadron. Its wave length is of order Aq CD , and 
so it cannot resolve the charge content of very much smaller bound states; in other 
words, it does not "see" the heavy quarks in a bound state of radius tq « Aq CD 
and hence cannot rearrange bonds. The resolution R(r) will approach unity for 
^Aqcd >> 1 and drop very rapidly with r for tAqcd < 1, in the functional form 

R(r) ~ (tAqcd)™, rA QC D<l, (4.29) 

with n = 2 [49] or 3 [39]. As a result, the integrand of eq. (4.28) will peak at 
some distance ro, with tq < r$ < Aqcd- Since the bound state radius of the 
quarkonium ground state decreases with increasing heavy quark mass, while R(r) 
is mQ-independent, ro — > tq — > as tuq — > oo. Hence P r vanishes in the limit 



rriQ — > oo because R(ro) does, indicating that the light quarks can no longer resolve 
the small heavy quark bound state. 



Fig. 4.3: Rearrangement transition p + J/ip — > D + D 



In a potential picture, the situation just described means that the charm quarks 
inside the J/ip have to tunnel from r = out to a distance at which the light quarks 
can resolve them, i.e., out to some r ~ cAq^ d , where c is a constant of order 
unity (Fig. 4.4). Such tunneling processes are therefore truly non-perturbative: 
they cover a large space-time region, of linear size ^qq D , and do not involve any 
hard interactions. Following [48], we shall here estimate the contribution of non- 
perturbative tunnelling to the dissociation of quarkonium states. 



Fig. 4.4: J/ip dissociation by tunneling 



In general, the problem of quark tunneling cannot be solved in a rigorous 
way, since it involves genuine non-perturbative QCD dynamics. However, the large 
mass of the heavy quark allows a very important simplification, the use of the 
quasiclassical approximation. In this approximation, the rate of tunneling -Rtun can 
be written down in a particularly transparent way: it is simply the product of the 
frequency ui^i of the heavy quark motion in the potential well and the tunnelling 
probability P t un when the quark hits the wall of the well, 

-Rtun = <^V-Ptun (4.30) 

The frequency ui^ is determined by the gap to the first radial excitation, 

lVtp ~ [M^i - Afy) ~ E^. (4.31) 

Consider now the potential seen by the cc (Fig. 4.4). For a particle of energy E, 
the probability of tunneling through the potential barrier V(r) is obtained from the 
squared wave function in the "forbidden" region. It can be expressed in terms of 
the action W calculated along the quasiclassical trajectory, 

Ptun = e~ 2W , (4.32) 

where 

W = \p\ dr, (4.33) 

J r\ 

Here the momentum \p\ is given by 

\p\ = [2M(V(x) - E)} 1/2 , (4.34) 

and ri, r2 are the turning points of the classical motion determined from the con- 
dition V{ri) = E. 

In our case, the width of the barrier is approximately 0.6 Aqq D , while its height 
(V — E) is equal to the dissociation threshold Eq. The mass M in Eq.(4.34) is the 
reduced mass, M = ijiq/2. We thus have 

W ~ 0.6 v/m Q E Q /A QCD . (4.35) 

For the J/ip, we get from (4.35) the value W ~ 3; this a posteriori justifies the use 
of quasiclassical approximation, which requires S > 1. 

Using eq. (4.35), we obtain as final form for the tunneling rate (4.30) 



.Rtun = E^ exp - (1.2y/m c E^/A QC d)- 



(4.36) 



With the above mentioned J/ip parameters, this leads to the very small dissociation 
rate 

R tun ~ 9.0 x 10" 3 fin -1 . (4.37) 
In terms of Ro, the J/ip survival probability is given by 

Stun = exp ^ - / dt R tun \ , (4.38) 



where t max denotes the maximum time the J/ip spends adjacent to the light hadron. 
In the limit t max — > oo, S^/ vanishes. However, the uncertainty relations prevent 
a localisation of the two systems in the same spatial area for long times. From 
Ax < Aqq D we get Ap > Aqcd, so that the longest time which the J/ip can spend 
in the interaction range of the light hadron is 



m 



= A QCD I 1 + T2— ) , (4-39) 
1 QCD 

with m for the mass of the light hadron. For nucleons or vector mesons, this 
time is 4 - 5 fm, and with this, the survival probability is very close to unity; 
hence non-perturbative tunneling interactions provide only negligible contributions 
to J/ip dissociation. 

In addition to such tunnelling, there can be direct and sequential thermal 
excitation to the continuum. In particular the latter still requires further analysis 
[48]. 

4.3 The Experimental Study of Charmonium-Hadron Interactions 

We now return to the break-up rate for J/ip interactions with slow light 
hadrons, as predicted by heavy quark QCD. Since this prediction is crucial for the 
use of quarkonia as confinement /deconfinement probe, it certainly must be checked 
experimentally. We shall now consider how that can be done. 

In Eq. (3.22), we had obtained the survival probability for charmonium states i 
in p — A interactions in terms of the break-up cross sections cij. We now just have to 
insert the cross section (4.24) into this expression to obtain the J/ip survival prob- 
ability as function of xf- In Fig. 4.5 the result is shown for a cms collision energy 
of 17.4 GeV, the value provided in CERN SPS P6-beam experiments. Included in 
this figure is also the colour octet interaction region xp > 0, together with data 
[27] and the fit obtained in section 3.2. We see that in the regime xp < —0.2, in 
which physical J/ip J s interact with the nuclear medium, the survival probability is 
essentially unity, due to the smallness of the cross section (4.24) in the threshold 
region. 



This picture has to be contrasted to the geometric absorption approach, which 
provides the basis for all hadronic accounts of .//^-suppression in nucleus-nucleus 
collisions. Here the J/ip cross section is assumed to attain its high energy value 
&j/iPn(s = oo) ~ 2.5 mb immediately at threshold. The survival probability for 
this case is also shown in Fig 4.5, and is seen to differ both qualitatively and 
quantitatively from that based on the QCD result (4.24). Through a measurement 
of J/^-production inp — A collisions at negative xp, technically attainable through 
data from a P6-beam incident on a hydrogen or deuterium target at positive xp, 
one can thus check if the threshold behaviour predicted by short distance QCD for 
inelastic J/^-hadron interactions is indeed correct. 

The actual experimental results will differ from what is shown in Fig. 4.5 for 
two reasons. We have already noted in section 3.4 the modifications expected 
because of the EMC effect. A second reason is that the J/ip peak observed in the 
measured dilepton spectrum is only to about 60 % due to directly produced 15" 
cc J/ip resonances; the remaining 40 % are mainly due to Xc production with the 
subsequent radiative decay Xc —> J/tp + 7 [22,23].* 

Concerning the EMC effect, we can either fold it into the predictions shown 
in Fig. 4.5, or it can be measured independently and then removed from the J/ip- 
production data. Measuring Drell-Yan dilepton production in the same kinematic 
region provides Ra/ p {x) (see Eq. (3.23)) for quarks directly, without any final state 
modification. A measurement of open charm production leads to Ra/ p (x) in just the 
same superposition of quark-antiquark annihilation and gluon fusion as in charmo- 
nium production, but again without any final state effect. Separate measurements 
of Drell-Yan and/or open charm production would thus determine the EMC mod- 
ification without addititional final state effects. Such measurements can therefore 
be used to unfold the EMC modification of J/ip data. 

To avoid any Xc admixture, it would of course be ideal to measure both J/ip and 
X c -production directly in p — A collisions [23] . This may be too difficult due to the 
abundance of photons in the case of heavy nuclei as targets. We shall therefore 
try to estimate the effect of the Xc admixture in the two scenarios considered here 
by simply adding the corresponding predictions with the noted 60/40 weights. For 
the geometric approach with its asymptotic cross sections, this requires a calcula- 
tion completely analogous to that for the J/ip, but now with a x ~ 6 mb and the 
correspondingly longer resonance formation time in the L of Eq. (3.22). 

The application of short distance QCD to calculate the inelastic x c -hadron 
cross section is as reliable as for the J/ip, since the binding energy of the Xc, £ x — 
0.24 GeV, is about equal to Aqcd- Nevertheless, such a calculation can give us 

* We shall for simplicity ignore here a further small contribution (< 5 %) fom 

ip' decay). 



some idea of the expected behaviour. The result is [25] 



2M x (m + e x ) 

(i - mi) 



6.5 



ct x at(s) ~ 11.3 mb x | 1 



(4.40) 



The asymptotic value is thus a factor two larger than the geometric estimate; this 
is a consequence of the fact that short distance QCD [40,41] leads to higher powers 
in the bound state radii than just r 2 . The behaviour of the \ C N cross section 
(4.40) is shown in Fig. 4.2. Inserting Eq. (4.40) into (3.22) and adding J/ip and 
Xc contributions then gives us the short distance QCD survival probability for the 
measured J/ip. It and the corresponding result from the geometric approach are 
shown in Fig. 4.6. As seen, for xp < —0.4, the two approaches differ qualitatively 
in their functional form and quantitatively by more than 20%. An experimental 
test of the theory for the interaction of heavy quarkonia with light hadrons should 
therefore be possible. 

As last point in this section, we comment briefly on the interaction of the 
ip' with light hadrons. Since its binding energy is only about 50 MeV, it lies almost 
at the open charm threshold and can definitely not be treated by short distance 
QCD. Here it might therefore not be unreasonable to assume that it attains its 
high energy value rather soon after threshold. This value can only be estimated by 
geometric arguments, and these suggest around 10 mb [46,47]. 

5. Quarkonium Dissociation in Confined and Deconfined Media 

Here we first show that at fixed temperature (or energy density) , a deconfined 
medium contains much harder gluons than a confined medium. Tightly bound 
quarkonia probe gluon hardness: while they were found to remain essentially un- 
affected in a confined medium at T < 0.5 GeV, a QGP of such temperature is 
shown to be very effective in their disssociation. To relate this to the environment 
produced in nuclear collisions, the resulting charmonium survival is studied in the 
case of isentropic longitudinal expansion. 

5.1 The Parton Structure of Confined vs. Deconfined Matter 

The ultimate constituents of matter are evidently always quarks and gluons. 
What we want to know is if these quarks and gluons are confined to hadrons or not. 
Let us therefore assume that we are given a macroscopic volume of static strongly 
interacting matter and have to analyse its confinement status. 

As prototype for matter in a confined state, we consider an ideal gas of pions. 
Their momentum distribution is thermal, i.e., for temperatures not too low it is 
given by exp(— E^/T) ~ exp(— p^/T). Hence the average momentum of a pion in 
this medium is (p^) = 3 T. The distribution of quarks and gluons within a pion is 



known from structure function studies; the gluon density is g(x) ~ 0.5(1 — x) 3 for 
large x = Pg/p^. 1 ^ As a consequence, the average momentum of a gluon in confined 
matter is given by 

1 3 

(Pg)conf = ^(Ptt) = ^T. (5.1) 

Hence in a medium of temperature T ~ 0.2 GeV, the average gluon momentum is 
around 0.1 GeV. In contrast, the distribution of gluons in a deconfined medium is 
directly thermal, i.e., exp(— p g /T), so that 

(p fl )deconf = 3T. (5.2) 

Hence the average momentum of a gluon in a deconfined medium is five times 
higher than in a confined medium 2 ); for T = 0.2 GeV, it becomes 0.6 GeV. An 
immediate consequence of deconfinement is thus a considerable hardening of the 
gluon momentum distribution. Although we have here presented the argument 
for massless pions as hadrons, it remains essentially unchanged for heavier mesons 
(p/co) or nucleons, where one can use a non-relativistic thermal distribution for 
temperatures up to about 0.5 GeV. We thus have to find a way to detect such a 
hardening of the gluon distribution in deconfined matter. 

The lowest charmonium state J/ip provides an ideal probe for this. It is very 
small, with a radius ~ 0.2 fm <C ^qcdi so that ^/V interactions with the 
conventional light quark hadrons probe the short distance features, the parton 
infra-structure, of the latter. It is very strongly bound, with a binding energy 
~ 0.65 GeV 3> Aqcd; hence it can be broken up only by hard partons. Since it 
shares no quarks or antiquarks with pions or nucleons, the dominant perturbative 
interaction for such a break-up is the exchange of a hard gluon, and this was the 
basis of the short distance QCD calculations presented in the previous chapter. 

We thus see qualitatively how a deconfinement test can be carried out. If we 
put a J/ip into matter of temperature T = 0.2 GeV, then 

- if the matter is confined, (p g ) C onf — 0.1 GeV, which is too soft to resolve the 
J/ip as a cc bound state and much less than the binding energy , so that the 
J/ip survives; 

- if the matter is deconfined, (p g )decon — 0.6 GeV, which (with some spread in 
the momentum distribution) is hard enough to resolve the J/ip and to break 
the binding, so that the J/ip will disappear. 

!) For very small x, recent results from HERA indicate a steeper increase; how- 
ever, this does not affect our considerations here. 

2 ) We could equally well assume matter at a fixed energy density, instead of 
temperature. This would lead to gluons which in case of deconfinement are approx- 
imately three times harder than for confinement. 



The latter part of our result is in accord with the mentioned prediction that the 
formation of a QGP should lead to J/ip suppression [1,16]. There it was argued 
that in a QGP, colour screening would prevent any resonance binding between the 
perturbatively produced c and c, allowing the heavy quarks to separate. At the 
hadronisation point of the medium, they would then be too far apart to bind to a 
J /t(j and would therefore form a D and a D. Although the details of such a picture 
agreed well with the observed J /if) suppression [50], it seemed possible to obtain 
a similar suppression by absorption in a purely hadronic medium [51], through 
collisions of the type 



Taking into account the partonic substructure of such hadronic break-up processes, 
we now see that this is in fact not possible for hadrons of reasonable thermal mo- 
mentum. Our picture thus not only provides a dynamical basis for J / ip suppression 
by colour screening, but it also indicates that in fact J/ip suppression in dense 
matter will occur if and only if there is deconfinement. We note, however, that 
the dynamical approach to J/ip suppression does not require a thermal equilibra- 
tion of the interacting gluons, so that it will remain applicable even in deconfined 
pre-equilibrium stages. 

While we have studied the hadronic part of the argument in detail in the previ- 
ous sections, we have so far not considered the dynamics of quarkonium dissociation 
in a deconfined medium. This will be taken up in the next section. 

5.2 Quarkonium Dissociation by Deconfined Gluons 

In section 4.2 we had obtained the cross section for the dissociation of a tightly 
bound quarkonium by an incident light hadron. Eq. (4.23) was obtained [2] by con- 
volution of the inelastic gluon-charmonium cross section with the gluon distribution 
in the light hadron. The gluon-quarkonium cross section itself is given by 



with k denoting the momentum of the gluon incident on a stationary quarkonium. 
In Fig. 5.1 we show the resulting break up cross section for gluon- J/ip and gluon- 
T interactions as function of the gluon momentum k. They are seen to be broadly 
peaked in the range 0.7 < k < 1.7 GeV for the J/ip, with a maximum value of about 
3 mb, and in the range 1.2 < k < 2.2 GeV for the T, with a maximum of about 
0.45 mb. Also shown in Fig. 5.1 are the corresponding cross sections for incident 
pions (note that now k = 3 in Eq. (4.2)), with high energy values of 3 mb and 0.5 
mb for J/ip and T, respectively. The break-up by incident hadrons is seen to be 
negligible up to momenta of around 4 GeV for the J/ip and 7 GeV for the T. Fig. 



J/ip + h^ D + D + X. 



(5.3) 




(5.4) 



5.1 thus provides the basis for the claim that in matter temperature T < 0.5 GeV, 
gluons of thermal momentum can break up charmonia, while hadrons cannot. We 
note here that, just as in the photoelectric dissociation of atoms, the break-up is 
most effective when the momentum of the gluon is somewhat above the binding 
energy. Gluons of lower momenta can neither resolve the constituents in the bound 
state nor raise them up to the continuum; on the other hand, those of much higher 
momenta do not see the (by their scales) large object and just pass through it. 

To illustrate this more explicitly, we calculate the break-up cross section for 
the J/ip as function of the temperature T of an ideal QGP. Using Eq. (5.4) with 
m c = 1.5 GeV and the J/ifj binding energy of 0.64 GeV, we then get 

a gJ/lp (T) = 65 mb x ^ . (5.5) 

The result is shown in Fig. 5.2 and confirms that up to about T ~ 0.5 GeV, only a 
deconfined medium can dissociate J/ip : s. We see moreover that the effective cross 
section for break-up in the temperature range 0.2 < T < 0.5 GeV is about 1.2 
mb. It is this value which will determine the suppression of the (pure 15) J/ip in 
a deconfined medium. 

Before we turn to charmonium production in a more realistic non-static envi- 
ronment, we want to consider briefly the possible role of quarks in the dissociation 
of charmonia in a QGP. This can be treated in a fashion quite similar to the interac- 
tion of quarkonia with light hadrons; the gluon distribution function in the hadron 
should now be replaced by an effective gluon distribution "in" quarks, i.e., the 
quark splitting function P(x) characterising the process q — > q + g. It's functional 
form for x — > 1 is fixed both by quark counting rules [52] and the Altarelli-Parisi 
equations [53], 

P(x)~(l-x) 2 , (5.6) 
which leads to an average gluon momentum 

3 

(^)deconf = 4 T - ( 5 - 7 ) 

The average momentum of gluons emitted by thermal quarks is thus higher than for 
gluons confined to thermal hadrons; it is nevertheless low enough to suggest that 
in an equilibrium plasma, the direct interaction with thermal gluons is the main 
dissociation mechanism. In pre-equilibrium, however, the quarks can be much faster 
and become the dominant cause of dissociation. 

5.3 Charmonium Survival in an Expanding Medium 

The exponential quarkonium survival probability of the general form (3.20) 
applies to a stationary medium of finite size. In nuclear collisions, the medium 



needs some proper time to to be formed, and it then expands until a time t / when 
the era of strong interactions in the medium ends. Hence slow quarkonia in the 
medium of not too high an initial temperature will in general stop interacting when 
the medium has cooled down enough, even though they have not yet left it. The 
survival probability of a quarkonium state % in such an expanding medium can be 
written as ^ 

S(T ) = exp{- [ 5 dt n{t) <n(t)}, (5.8) 

J t 

where n(t) is the density of scattering centers at time to < t < tf and (Ti(t) the 
break-up cross section for state % in the medium at that time. If we assume isentropic 
longitudinal expansion, density, temperature and time are related by 

n(t ) _ (To) 3 = l 
n(t) \TJ t ' 

with T < Tq. Using this relation, we can rewrite Eq. (5.6) in the form 

f Tc a h (T) f T ° a 9 (T) 

S(T ) = exp{-3n h (T )t / dT^-1 - 3n g (T )t Q dT^-±}, (5.10) 

JT f 1 JT C 1 

where Tf is the temperature of the medium for which the break-up of state i stops. 
The first term in the exponent corresponds to dissociation by hadrons, the second 
by gluons. We had seen that the main contribution to the dissociation comes from 
gluon-quarkonium interactions, making S(T ) ~ 1 up to T ~ T c ~ 0.15 GeV. 
Hence by inserting the cross section (5.5) (Fig. (5.2)) into Eq. (5.10), we obtain the 
survival probability of a J/ip in a medium in isentropic longitudinal expansion, for 
an initial temperature Tq > T c ; below that, it is unity. The thermalisation time to 
is generally argued to be around 1 fm; and for an ideal QGP, the density of gluons 
at temperature T is given by the Stefan-Boltzmann form 

n(T) = ^T 3 . (5.11) 

The resulting suppression as function of temperature is shown in Fig. 5.3. The 
survival probability is seen to decrease very rapidly with increasing temperature, 
essentially vanishing for T > 0.3 GeV. 

6. Pre-Equilibrium Deconfinement 

In the previous chapter, we established quarkonium suppression as a probe 
to check if a given sample of strongly interacting matter consists of deconfined 
quarks and gluons. In an actual nuclear collision, however, such a suppression 



(5.9) 



could have been caused before the onset of thermalisation. Here we therefore study 
the possibility of pre-equilibrium quarkonium suppression. 

6.1 Shadowing and J/ip Suppression in Nucleus-Nucleus Collisions 

Consider a nucleus-nucleus collision at CERN SPS energy (y/s ~ 20 GeV in 
the nucleon-nucleon cms), leading to J/ip production at mid-rapidity, yj/^=0, with 
the production mechanism as described in Chapter 2. During the first 0.25 fm, 
the produced cc pair is a colour octet and will interact as such with the passing 
nucleons of target and projectile. The momenta of these nucleons in the rest system 
of the J/ip will be around 10 GeV. The situtation seen by the J/ip is thus very 
similar to that encountered at xp = in collisions of 200 GeV/c protons with a 
nuclear target. As discussed in Chapter 3, one here finds a suppression dominated 
by "nuclear shadowing", i.e., destructive interference of the scattering on different 
nucleons. Such an effect will now arise from both projectile and target, however, 
and this suppression must be removed before final state effects can be studied [11]. 

To see the effect of this shadowing correction, we show in Fig. 6.1 the latest 
data taken by the NA38 collaboration at CERN in p — U and S — U collisions 
[54,55]. Here the suppression is measured with respect to the high mass Drell- 
Yan continuum, isospin corrected for p — U. The nucleus-nucleus data are shown 
as function of the neutral transverse hadronic energy Ej, produced in the collision. 
The p—U value is simply included in this figure; it is not associated to any particular 
Ej,. The difference between the p — p and p — A results is, as discussed in Chapter 
3, essentially given by the shadowing function R s ^/ p {xf — 0). To correct the data 
for shadowing effects, we therefore divide the p — U data by Rfjy p {0) and the S — U 
data by R s s h / p {0) x R^ p {0). The result is shown in Fig. 6.2; the p — U value and that 
for the S beam at low Ej, now agree. The remaining suppression in the nuclear 
collision data is now due to effects on the colour singlet cc in its state after the time 
of colour neutralisation. 

We thus want to study the effect of the environment on this cc. If we ignore 
nuclear stopping, target and projectile nucleons retain in the cc cms their initial 
momenta Pn — 10 GeV/c. The cc then sees two nuclei, each Lorentz-contracted to 
about 1 fm, since {2m/ y/s) ~ 0.1; for simplicity, we consider here A — A collisions. 
At the colour neutralisation time rg ~ 0.25 fm, the centers of these thin discs 
have become separated by approximately 0.5 fm, so that the target and projectile 
nucleons still have considerable overlap. Stopping would slow down the nucleons 
to increase this. The detailed kinematics encountered by the colour singlet cc is, 
however, quite unimportant for the essential question: can a pre-equilibrium state 
of hadrons account for the observed suppression? The initial medium, at the point 
of maximum overlap between target and projectile, is one of twice nuclear density, 
containing nucleons of 10 GeV/c momentum in the cc rest system; equivalently, we 



can assume normal nuclear density and an effective cc path length L = 2 x 3Ra/4=, 
with Ra = 1.12 A 1 / 3 [34]. Whatever scattering processes now occur lead to an 
increase of the density n, but at the same time to a decrease of the momenta P of 
the hadronic constituents of the medium. Momentum conservation requires for the 
momentum flow through a given surface F 

P n F = P N n F = const. (6.1) 

so that the momentum in fact drops as 1/n. If we neglect the cross section reduction 
during the time needed for the singlet cc to become a full-sized resonance, the 
survival probability of a J/ip in the medium becomes 

Sj/V = exp{-n a(P N )L}. (6.2) 

The medium actually seen by the J/ip will have a higher density [56], but its con- 
stituents will have lower momenta. In view of Eq. (6.1), we can nevertheless use the 
form (6.2) to calculate the survival probability. In contrast to the phenomenolog- 
ical fit of [34], Eq. (6.2) contains the break-up cross section at the actual collision 
energy, and in Chapter 4 we have seen that this cross section is much below its 
geometric high energy value. At Pjv = 10 GeV/c, we find (see Fig. 4.2) Gj/^ ~ 0.8 
mb. Inserting this into Eq. (6.2), with an average value L ~ 10 fm for central 
S — U collisions, leads to S ~ 0.9 for the survival probability. Since we have here 
neglected the cross section reduction during the J/ip evolution as well as the path 
length reduction by the part already included in the shadowing correction, the ac- 
tual survival probability will be larger. Thus hadronic pre-equilibrium interactions 
cannot account for the measured J/ip suppression. 

We therefore conclude that the J/ip suppression observed in the NA38 exper- 
iment [15] provides evidence for the existence of deconfined gluons in the medium 
probed by the J/ip. By this we mean that the medium in which the J/ip finds itself 
after a central nucleus- nucleus collision contain gluons whose momentum distribu- 
tion is harder than that found in mesons or nucleons [57]. Only in collisions with 
such gluons can the tightly bound J/ip be dissociated; gluons confined to hadrons 
are not sufficiently hard. We emphasize that this conclusion makes crucial use of 
the energy dependence of the J / flight hadron break-up cross section as calculated 
in short distance QCD. As pointed out in Chapter 4, this result can and should be 
experimentally checked. 

6.2 Deconfined Gluons and QGP 

The existence of deconfined gluons is clearly not equivalent to the existence 
of a QGP, in which such gluons are in thermal equilibrium. It is only the first 
step towards a QGP: it shows that in the large-scale medium there exist deconfined 



partons. Their thermalisation comes at the end of a parton interaction cascade, 
and it is not at all clear whether at present energies there are enough deconfined 
partons in a sufficiently large and long-lived system to reach this stage. How can 
we check experimentally whether parton thermalisation is or is not achieved? 

The cascade formation of a QGP [58,59] starts with the production of gluons 
in primary collisions; these then interact again to produce secondary gluons, quite 
possibly through multigluon production [60], and so on, until production and ab- 
sorption balance to form an equilibrium system. In equilibrium, the number of 
gluons per unit volume is simply determined by the energy (or entropy) density 
of the system. In the pre-equilibrium stage, it is lower and (still) proportional to 
the number of primary collisions. We thus expect J/ip suppression by deconfined 
gluons in a pre-equilibrium system to increase with the number of nucleon-nucleon 
collisions (or equivalently, with the effective path length of the J/ ip) . On the other 
hand, equilibrium suppression would be independent of the number of such colli- 
sions and depend on the effective energy density of the system only [57]. 

To test this, we can study Pb — Pb collisions as function of increasing cen- 
trality (i.e., of increasing Et)- In this case, the energy density remains essentially 
unchanged [29] while the number of collisions increases.* If the J/ip suppression is 
found to increase, the system cannot have reached equilibrium, and hence the sup- 
pression must be due to deconfined gluons in the pre-equilibrium stage. Once the 
suppression becomes independent of the number of primary collisions, equilibrium 
and hence the QGP is reached. The forthcoming NA50 P6-beam data from the 
CERN-SPS should thus be able to resolve this question. 

7. Conclusions 

We first summarize the essential theoretical points of this work, and then survey 
the main conclusions we think can be drawn from present data. 

Quarkonium production in hadron-hadron collisions is today quite well under- 
stood in terms of elementary parton interactions (gluon fusion, quark-antiquark 
annihilation). The distributions of the partons within the colliding hadrons are 
determined in deep inelastic lepton-hadron scattering, and the (non-perturbative) 
colour neutralisation of the produced heavy QQ pair can also be fixed empirically. 

Hadron-nucleus collisions determine what happens to quarkonium production 
in a confined medium. In the production of fast quarkonia in the nuclear rest frame, 
a fast colour octet QQ passes the medium, leading to quantum-mechanical inter- 
ference ( "nuclear shadowing" ) and energy loss. Slow quarkonia in the nuclear rest 

* Since a change in centrality of S — U collisions changes both the number of 
primary collisions and the effective energy density, this does not distinguish pre- 
equilibrium from equilibrium. 



frame are subject to EMC effect modifications of the colliding partons in addition 
to collisions with nucleons in the nucleus. 

Because of the small size and the large binding energy of the lowest quarkonium 
states, their interaction with light hadrons is calculable in short distance QCD. 
They can interact in leading order only through the exchange of a hard gluon, 
and the gluon distribution in the light hadrons is known to be very soft. The 
resulting prediction is a cross section which rises very slowly from threshold to its 
high energy value, suppressing strongly any break-up of quarkonium ground states 
by slow mesons or nucleons. 

As a consequence of this suppressed quarkonium dissociation by light hadrons, 
confined matter (at temperatures T ~ 0.5 GeV) becomes transparent to J/ip's. 
The momentum of deconfined thermal gluons, on the other hand, is large enough to 
give rise to effective J/ip dissociation; such dissociation can occur also by deconfined 
gluons not in equilibrium. Strongly interacting matter thus leads to J/ tf) suppression 
if and only if it is deconfined. The loosely bound ip' can be broken up in both 
confined and deconfined matter, though presumably more in a deconfined medium. 

In nucleus- nucleus collisions, the interaction of the coloured cc pair with target 
and projectile nucleons can be taken into account through the nuclear shadowing 
determined in hadron- nucleus collisions. Any strong J/ip suppression remaining 
after removal of these shadowing effects is accountable only by interaction with 
deconfined gluons. 

What have we then learned from the h — A [27,28] and A — B [15,61] data 
available so far? What should the forthcoming P6-beam data clarify? 

The equality of J/ip and ip' suppression in h — A collisions for xf > 0, as 
well as the size and xf dependence of the observed effect, are in full accord with 
the passage of a colour octet through nuclear matter. The equality of J/ip and 
ifj' suppression and the observed xf dependence are in clear disagreement with 
any description based on the absorption of physical charmonium states in nuclear 
matter. 

There is a lack of data for charmonium production in a kinematic regime in 
which fully formed J/^'s could interact with nuclear matter. Such data could be 
obtained by experiments using the P6-beam incident on a light target [25]. 

The observed difference between J/tjj and ip' suppression in nucleus-nucleus 
interactions [61] indicates that the relevant mechanism here is different from that 
in h — A collisions. 

The observed J /if) suppression cannot be accounted for in terms of nuclear 
shadowing on both projectile and target; even after removal of shadowing, an in- 
dependent suppression (about 40 % between low and high Et) remains. 

If the strong threshold suppression of J/ip break-up by light hadrons, as pre- 
dicted in short distance QCD, is experimentally confirmed, the J/tjj suppression 



observed in O — U and S — U interactions can only be accounted for by the pres- 
ence of deconfined gluons. 

It remains open whether the deconfined gluons required for J / ip suppression in 
nuclear collisions are already equilibrated; hence their existence does not establish 
QGP formation, but quite likely only a first step towards a thermalised deconfined 
medium. If the gluons are not in equilibrium, the resulting J/ip suppression should 
increase with increasing Et in Pb — Pb collisions at fixed energy; in equilibrium, 
the suppression would remain approximately constant. 
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